Hybrid multi-site excitations in dipolar condensates in optical lattices 
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Strong ID lattices usually lead to unconnected two-dimensional gases. The long-range character 
of the dipole-dipole interactions leads to a novel scenario where non-overlapping gases at different 
sites may interact significantly. We show that the excitations of non-overlapping condensates in 
fD optical lattices acquire a band-like character, being collectively shared by different sites. In 
particular, the hybridization of the modes significantly enhances the rotonization of the excitations, 
and may induce roton-instability. We discuss the observability of this effect in on-going experiments. 
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A novel path on cold gases is currently being opened 
by recent experiments in which (magnetic or electric) 
dipole-dipole interaction (DDI) plays a significant or even 
dominant role. On one side the recent creation of het- 
eronuclear molecules in the lowest ro-vibrational level 
P, 0] opens exciting perspectives for the achievement 
of a quantum degenerate gas of polar molecules, which 
may possess large dipole moments (e.g. ~ 0.5 Debyes 
for KRb [H). On the other side, the magnetic DDI has 
been shown to lead to exciting novel phenomena in re- 
cent experiments in Bose-Einstein condensates (BECs) 
of Chromium (which has a magnetic moment fi — 6/iB, 
with /is the Bohr magneton) [3l|. Particularly interesting 
is the fact that the short-range interactions (SRI) may be 
suppressed by means of Feshbach resonances, leading to 
a purely dipolar gas Q. The DDI plays also a signifi- 
cant role in very recent experiments on spinor Rubidium 
BECs, in spite of the small /i — l/^s, since the energy 
scale of the DDI becomes comparable with the (also very 
small) energy scale of spin-changing collisions Very 
recent experiments have shown as well that the DDI leads 
to a observable damping of Bloch oscillations in Potas- 
sium BECs in tilted optical lattices [6]. 

The partially attractive character of the DDI results 
in nontrivial stability conditions for dipolar BECs. Low- 
momentum instability (phonon instability) results in 
a geometry-dependent instability against collapse in 3D 
traps, as recently observed experimentally Q, or soliton 
formation in 2D geometries '9l|. Interestingly, the mo- 
mentum dependence of the DDI allows for a second type 
of instability (roton instability) related to the appear- 
ance of a roton-like minimum in the dispersion law of 
elementary excitations |1Q] . Roton instability leads to 
local collapses [11] or stabihzed modulated density pro- 
files in sufficiently tight traps [l^ . 

The long-range character of the DDI induces a nonlo- 
cal nonlinearity in dipolar BECs that resembles that en- 
countered in plasmas [l^] or nematic liquid crystals [l^ . 
As a consequence, novel phenomena as stable 2D soli- 
tons become possible [H, ll^]- This nonlocality leads to 
fundamentally new physics for quantum gases in optical 
lattices, since it induces interactions between neighboring 



Figure 1: Scheme of the system under consideration. 



sites. As a consequence, dipolar bosons in optical lattices 
are described by extended versions of the Bose-Hubbard 
Hamiltonian, and may present a wealth of novel phases, 
as supersolid [13] or Haldane-phases [11]. In addition, 
contrary to the case of SRI, very deep optical lattices 
do not lead to independent low-dimensional gases, since 
non-overlapping atoms at different sites interact. As a 
consequence nonoverlapping BECs in two- well potentials 
may scatter Il9l . pair superfluidity may appear in ladder- 
like lattices [20[, and even filament condensation may 
occur [2l|. The effects of the intersite DDI have been 



observed experimentally for the first time in very recent 
experiments in Florence on Bloch oscillations 

This Letter is devoted to the analysis of non- 
overlapping dipolar BECs placed at different sites of a 
deep two-well potential or ID optical lattice. As men- 
tioned above, contrary to the case of purely SRI, the 
deep potential does not lead to independent 2D BECs. 
In particular, we show that the elementary Bogoliubov 
excitations of disconnected BECs placed in a two- well po- 
tential couple through the DDI leading to hybrid modes 
which are collectively shared by both wells. Interestingly 
this hybridization may significantly alter the stability of 
the system against roton instability. We show that this 
effect is significantly enhanced for the case of a ID optical 
lattice with multiple sites, where a band-like spectrum is 
induced by the inter-site DDI. We analyze in detail the 
stability diagram, and finish with a discussion of the ex- 
perimental observability in different experiments. 

In the following, we consider a BEC of particles with 
mass m and electric dipole d (the results are equally 
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valid for magnetic dipoles) oriented along z by an ex- 
ternal field, and that hence interact via a dipole-dipole 
potential: Vrf(r) — (/^(l — 3cos^(6'))/r'^, where B is the 
angle formed by r with the z axis. At sufficiently low 
temperatures the physics of the dipolar BEC is provided 
by a non-local non-linear Schrodinger equation (NLSE): 



-^V2 + ^(f)+5|*(r,t)P 



(1) 



where g = iirH'^a/m, with a the s-wave scattering length, 
and V{r) is the trap potential. In the following we use 
the convenient dimensionless parameter /3 = gd/g, that 
characterizes the strength of DDI compared to the short 
range interaction, where gd — STrd^/S. Note that P may 
be easily controlled experimentally by means of Feshbach 
resonances, as recently shown in Ref. [§]. 

We consider first the case of a quasi-2D homogeneous 
BEC confined in z by V{r) = muj'^z'^ /2, which is suffi- 
ciently strong such that ^'(r, = ±{p,t)^o{z), where 
$0(2) = exp(-z^/2^^)/7r^/^/ y^ is the ground state of the 
transversal oscillator (l^ = y^h/muiz)- The ground state 
of the homogeneous 2D BEC is of the form ^'^(p, i) = 
exp{—i{^/h + LUz)t)^/no, where no is the 2D density, and 
fj, is the 2D chemical potential. Introducing this form 
into the NLSE ([T]), one obtains /i = {g + gd)no/V2TTlz 
(the 2D condition is satisfied for fi ^ ^z)- Inserting 
a plane-wave Ansatz 4'(r, t) = {y^no + Wgexp(i(f • p — 
iet/h) — V* exp{~iq ■ p~ iet / h))^o{z) exp{—i{fi/h + uJz)t) 
into the NLSE ([T]), and linearizing in Uq, Vq, we obtain 
the Bogoliubov spectrum of elementary excitations: 



e{q) = {Eq [Eq + 2A]} 



1/2 



(2) 



where Eq = h^q^/2m, A = p. - {gdno/V2nlz)F{qlz/V2) 

and F{x) = ^|a;|erfc(x)e"''. Note that without DDI 
{f3 — 0) we recover the usual Bogoliubov spectrum for 
a 2D BEC with purely SRI. In particular, if a < and 
P — 0, e{q)'^ < for g — !■ 0, recovering the well known 
phonon instability (and subsequent collapse) in homoge- 
neous BEC with a < 0. If the dipolc is sufficiently large, 
such that g + gd > 0, then the DDI prevents the insta- 
bility at g ^ 0. However, due to the g-dependence of 
the DDI (given by the monotonously decreasing charac- 
ter of the function F), the dispersion e{q) may show for 
intermediate gd values a roton-like minimum at a finite 
value ofqlz (Fig.[2]). For sufficiently low DDI e{q)^ < at 
the roton-like minimum, leading to dynamical instability 
(roton instability). For > Per (with Per dependent on 
the ratio gno/lzfiLLiz) roton instability is prevented, and 
the 2D homogeneous BEC is stable. 

In the following we show that Per is significantly modi- 
fied in the presence of other neighboring quasi-2D dipolar 
BECs. We consider the case of an optical lattice along z 
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Figure 2: Roton minimum in the dispersion relation e/hcOz of 
a single 2D BEC with /3 = -1.07, L = 0.09pm, a = -2 nm, 
and a 3D density of no/V^I-z = lO^'^cm"^. 

(Fig.[T]) described by a potential I^(r) = sEr sin^(7rz/A), 
where A is the intersite spacing, and s provides the lat- 
tice depth in units of the recoil energy Er = fi.^7r^/2mA^. 
As in the previous discussion we consider no trapping on 
the xy-plane (we discuss the potentially important role 
of the harmonic confinement on the xy-plane at the end 
of this Letter). At each lattice node V{r) may be ap- 
proximated by an effective harmonic oscillator potential 
Veff{z), with effective oscillator length Iz « As~^^'^/n. 
The lattice is considered strong enough so that we can 
assume that there is no spatial overlap between wave- 
functions in different lattice sites, and hence we may ne- 
glect any hopping. We assume the DDI small enough to 
neglect pairing [2o'| or filamentation [2l| . 

We start our discussion on inter-site effects with the 
two- well case. This simplified scenario already captures 
many features of the effect discussed. In addition, two- 
well potentials may be experimentally realized and are 
currently of considerable interest 2^, l23| . The quasi-2D 
BEC in the i-th layer is given by the extended NLSE: 



iS^*,(r,t) = 



-^V' + Veffiz)+g\^,{r,t)\^ 



*(r,0, 



(3) 



where for a two-well potential, i,j = 1,2. Note that, 
crucially, the DDI couples now the i-th layer to the j- 
th one. Similar to the single-site discussion, we consider 
a strong z-confinement at each site, and hence we may 
employ a quasi-2D Ansatz 5'i(r) — ^'^,i(/3, i)<f>o(-z — Zi), 
where $o(^) has the form discussed above, and Zi is 
the position of the i-th lattice node. The ground-state 
of the condensates at the two layers in given by the 
Ansatz 5'j_,i(p, i) = y^e~*^^/''+"')*, where we consider 
the same 2D density no at both sites. Introducing this 
Ansatz into the NLSE ([3]) we obtain the 2D chemical po- 
tential p = /i-|-A(A), with /i the chemical potential of an 
individual well and A(A) = {gdno/V2Trlz)e~^ /'^ 

As above we are interested in the elementary excita- 
tion of these systems. For A ^ 00 the Bogoliubov modes 
at each site are independent and described by the single- 
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site expression ([2]). For finite A the inter-site coupling 
leads to an hybridization of the modes at both sites with 
significant consequences, as we discuss below. As for 
the single-site discussion we insert a plane-wave Ansatz 
*i(f,i) = {y^ + Uq^exp{iq- p-iet/h)-v*^exp{-iq- p~ 
i€t/h))^o{z) exp{-i{jl/h + uj;,)t) into the NLSE and 
linearize in Uqi, Vqi. In this way we obtain four coupled 
Bogoliubov-de Gennes equations for 2,^1,2}, which 
may be diagonalized to obtain the Bogoliubov modes: 

e±{q) = {Eq[Eq + 2A±CiA)]}'/\ (4) 



where 




Figure 3: Dispersion law (in units of £"0 ~ h/mA^) for a single 
site (dashed) and for two wells (solid) for /3 = —1.2, A = 0.53 
fj,m, s = 13.3, a = —2 nm, and no/\/2Trlz — 10^''/cm^. 

Note that for A — ^ 00, C(A) — and we recover two 
degenerate independent modes. For finite A the modes 
at the two wells hybridize, and two different branches ap- 
pear for each q, one stiffer than the modes for A 00, 
and the other softer. The latter is particularly interest- 
ing, since the soft mode is more prone to rotonization 
(Fig. [3]). Interestingly, under proper conditions, two par- 
allel non-overlapping BECs may become roton-unstable 
even if they were stable separately. As a consequence, a 
larger Per is necessary to stabilize the two-well system. 

The hybridization (and consequent destabilization) in 
two-well potentials becomes even more pronounced for 
the case of dipolar BECs at A^s > 2 sites of a ID optical 
lattice, since a site i couples with all its neighbors j (of 
course with decreasing strength for growing |i — j|). For 
simplicity of our analysis we consider the case in which 
all lattice sites present the same 2D density Uq. In that 
case, one may easily generalize the two-site analysis to 
the multi-site case, to reach a set of coupled Bogoliubov- 
de Gennes equation for fqi = Uqi + Vqi : 

e^fq, = Eq{Eq + 2A)fq, + 2Eq ^ C(A|z - j\)fqj. (6) 
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Figure 4; Band-like dispersion (in units of Eo = ft/mA^) for 
Ns — 40 and /3 = —2.44. Other parameters are as in Fig. |3] 
We indicate the dispersion law for Ng — 1. 

After diagonalizing the matrix of coefficients at the rhs of 
Eqs. (Ill), we obtain numerically the corresponding band- 
like set of Ns elementary excitations (Fig. |4]). Note that 
the band-like spectrum has an upper phonon-like bound- 
ary which for large has an approximate sound velocity 
Cs ~ ^(A-|-X;„ C'(A|n|))/TO. The lower mode of the N, 
manifold becomes significantly softer than the individual 
modes for independent sites. As a consequence the ro- 
ton instability extends to larger (3cr when Ng increases, 
until saturating for a sufficiently large Ns (due to the 
decreasing DDI for increasing distance between sites). 
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Figure 5: Stable and unstable regimes for Ns dipolar 2D 
BECs. We employ the same parameters as in Fig. [S] 

Fig. [S] summarizes our results on the stability as a func- 
tion of f3 (wc recall that g < 0). As mentioned above 
ifff + ffd < (|/3| < 1) the system is unstable against 
phonon instabiHty. For 1 < < \f3cr{Ns)\ the system 
is unstable against roton instability. \Pcr\ increases when 
Ns grows until saturating for sufficiently large Ns- For 
\I3\ > \l3criNs)\ the quasi-2D BECs are stable. 

The value of qrot when the roton becomes unstable is 
of particular importance. Fig. [6] shows a typical varia- 
tion of qrot A at the curve /3 = ficr{Ns) as a function of 
Ns- Note that qrot at Per shows a maximum for small 
Ns- For small Ng, Per (and hence the on-site repulsive 
DDI) increases significantly, and hence the value of qrot 
at Per increases. For larger Ns Per tends to saturate, as 
mentioned above, and the repulsive on-site DDI remains 
approximately constant along the curve Pcr{Ns). As a 
consequence, the increase in Ns just increases the attrac- 
tive contribution of the DDI of neighboring sites, and the 
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DDI becomes less effective in compensating the attrac- 
tive on-site SRI. As a result of that, qrot decreases until 
saturating at a value lower than that for a single site. 
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Figure 6: Roton momentum qrot when the roton minimum 
touches zero as a function of the number of sites Ns- We 
employ the same parameters as in Fig. O 

Typical experiments work with an harmonic xy- 
trapping (of frequency oJxy)- Although we have assumed 
homogeneous quasi-2D gases, we may estimate the ef- 
fect of the xy-trapping by considering an effective cut-off 
at low momenta q^ut — ^/Ixy, where l^y — \J h/niujxy 
is the harmonic oscillator length characterizing the xy- 
trap. In a good approximation we may consider that all 
features occurring at momenta q < qcut are suppressed 
by the trap. As a consequence one expects that the xy 
confinement suppresses roton instability for frequencies 
ujxy > i^cut- For typical densities lO^^cm"^, and typical 
intersite separation A = 0.53/im, we estimate for ^^Cr 
that for a single site Per is achieved at a scattering length 
a ~ — 31ao, and that for this case LOcut — 66Hz. For the 
same case but Ng = 4 (which is the maximum of the 
corresponding g^ot curve), (3cr is achieved for a ~ — 24ao, 
and oJcut — 160Hz. For the latter case an instability rate 
of ~ 5ms is expected for a = — 24.5ao. For ^^K, the 
numbers are more restrictive (due to the lower magnetic 
moment). For Ng = 25 (maximum of the qrot curve), f3cr 
is achieved for a ~ — 0.52ao, and Ucut — 4Hz. For this 
case, one expects F"^ ~ 180ms at a = — 0.53ao. 

Note that the fact that qrot shows a maximum may 
have interesting consequences in experiments, since this 
suggest that for some intermediate trapping frequencies 
the instability may be just present for a given window of 
values of Ng. Note also that the previous discussion just 
refers to the destabilization when the roton touches zero. 
For even larger values of |a|, a larger region of q may be- 
come unstable. However if the xy-trap just allows for the 
resolution of the upper boundary of the unstable region 
and not the lower one, roton and phonon instability may 
become experimentally undistinguishable. 

Summarizing, the nonlocal character of the DDI leads 
to a novel scenario where non-overlapping gases at dif- 
ferent sites interact significantly. Contrary to the case of 
pure SRI, the DDI leads to the hybridization of the exci- 
tations at different sites, which acquire a collective band- 
like character. In particular, the hybridization of the 



modes leads to a significant enhancement of the rotoniza- 
tion of the excitations, and may induce roton-instability 
for values of the SRI at which a single site is stable. Fi- 
nally, we have discussed the experimental requirements 
for the observation of the roton instability. 

Note: After the completion of this work we became 
aware of a similar analysis by Wang and Demler 
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which roton-softening due to intersite interactions is dis- 
cussed in the context of recent experiments in Florence 
0. Although we consider that roton softening plays no 
significant role in the damping observed in Ref. 6] due 
to the a;?;-confinement, a weaker u^y (along the lines dis- 
cussed above) could allow for the instability discussed by 
Wang and Demler, and by us in this manuscript. How- 
ever, a more careful quantitative analysis is necessary, 
taking into account both the xy-trapping, and the z- 
trapping, which we plan to investigate in a further work. 
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